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Abstract. - We report the observation of dynamical behaviour in dilute, aqueous solutions of 
a surfactant CTAT (cetyl trimethylammonium p-toluenesulphonate) , below the overlap concen- 
tration c*. At these concentrations, CTAT forms cylindrical micelles and shows a pronounced 
shear thickening transition above a concentration-dependent critical shear rate j c . An analysis 
of the time-series of the stress relaxations at controlled shear rates in the shear-thickening regime 
shows the existence of correlation dimensions greater than two and positive Lyapunov exponents. 
This indicates the existence of deterministic chaos in the dynamics of stress relaxation at these 
concentrations and shear rates. The observed chaotic behaviour may be attributed to the 
stick-slip between the shear - induced structure (SIS) formed in the sheared surfactant solution 
and the coexisting dilute phase. At still higher shear rates, when the SIS spans the Couette, there 
is a transition to higher-dimensional dynamics arising out of the breakage and recombination 
of the SIS. 



Introduction. - It is well known that surfactant molecules, under appropriate conditions, 
can assemble reversibly in solution to form various supramolecular structures. In very dilute 
solutions, at concentrations above the critical micelle concentration, these molecules self- 
assemble to form spherical micelles. Increasing the surfactant concentration results in the 
formation of cylindrical structures, which overlap and entangle to form clear viscoelastic 
solutions of giant wormlike micelles. At still higher concentrations, the surfactant solutions 
form mesophases exhibiting nematic, hexagonal or lamellar ordering. It is possible to pass from 
one phase to the other by simply changing the surfactant and cosurfactant concentrations, the 
salinity or the temperature of the system. The cylindrical structures that are formed at very 
low concentrations c (c < c* , c* is the overlap concentration), get entangled on shearing at 
a shear rate 7 > j c to exhibit a pronounced increase in viscosity jl], [| [|. At 7 < 7 C , the 
surfactant solutions exhibit Newtonian flow behaviour. Recently, direct imaging studies of light 
scattered by the shear-thickening CTAB/ NaSal and TTAA/ NaSal have shown the existence 
of unstable shear induced structures/ phases (SIS/ SIP) which fracture at very high 

shear rates, giving rise to shear thinning H. Direct imaging of the SIS using freeze- fracture 
electron microscopy shows that these structures span hundreds of microns and have sponge- like 
textures [Q. 
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In this paper, we present an experimental study on the nonlinear dynamics of stress 
relaxation in dilute solutions of the surfactant CTAT, at concentrations below the overlap 
concentration c*. At these concentrations, CTAT forms cylindrical micelles || and undergoes 
a shear-thickening transition when sheared at 7 > %, which is followed by shear-thinning 
at higher values of 7 ED]. The phase diagram and flow behaviour of aqueous solutions of 
CTAT in the linear and nonlinear regimes have been investigated extensively by Soltero et al. 
[0. Previous investigations of the flow behaviour of shear-thickening solutions of CTAT [[[] 
using rheometry and neutron scattering have shown the coexistence of a highly viscoelastic 
shear-induced phase (SIP) with a viscous regime made up of short aggregates. Oscillations 
in the apparent viscosity (and hence the stress) of dilute CTAT at controlled shear rates in 
the shear thickened regime have been observed by Gamez-Corrales et al. [0. Observations of 
viscosity oscillations in CTAB/NaSal and TTAA/NaSal H in water under controlled stress 
conditions have been explained in terms of the growth and retraction of the SIS that form in 
the sheared solution. 

In a recent paper jfjl, we have established the existence of deterministic chaos |?J in the 
stress relaxation of shear-thinning CTAT solutions at concentrations c > c*. The oscillations 
in the normal and viscoelastic stresses were explained in terms of the stick-slip between shear 
bands Q that form in the plateau region of the flow curve. Here, we extend our experiments to 
more dilute solutions of CTAT (c = llmM, zero-shear viscosity rj ~ 2mPas, c* ~ 15mM) and 
show that the occurence of coexisting phases in a sheared surfactant solution is a prerequisite 
for the existence of dynamical behaviour in its stress relaxation. In the following sections, we 
will demonstrate the existence of deterministic chaos in the stress relaxation of shear-thickened 
CTAT, possibly due to the stick-slip between the SIS with the coexisting dilute phases. To our 
knowledge, this is the first observation of deterministic dynamics in the stress relaxation of 
shear-thickening surfactant systems. The observed chaotic behaviour is followed by an increase 
in complexity of the dynamics of stress relaxation at higher shear rates due to the fracture of 
the SIS, which is accompanied by the formation of large vortices in the sheared solution 

Sample Preparation and Experimental Apparatus. - CTAT, purchased from Sigma Chem- 
icals, Bangalore, India was used to prepare the samples by dissolving requisite amounts of 
CTAT in doubly distilled and deionised water. The sample containers were sealed and kept 
at 25° C and shaken frequently to accelerate homogenisation. By measuring the viscosities of 
CTAT solutions prepared at different concentrations as a function of the shear rate 7, and 
extrapolating the values to 7=0, we have estimated the overlap concentration c* of CTAT in 
water to be ~ 15mM. The rheological measurements were made using a Rheolyst AR1000N 
stress-controlled rheometer at T=25°C, in a concentric cylinder geometry with an conical end 
(outer cylinder diameter = 30mm, inner cylinder diameter = 27.66mm, immersed height = 
25.8mm). 

Experimental Results. - 

Figure 1 shows the flow curve of llmM (0.5wt.%) CTAT at T=25°C, measured by allowing 
30 seconds between the acquisition of successive data points. The flow curve shows all the 
distinct regimes seen in TTAA/NaSal by Hu et al. in stress-controlled experiments ft], viz., 
a Newtonian region of constant viscosity (Regime I) followed by a shear thickening regime 
where the SIS nucleates heterogeneously from the inner wall of the Couette cell (Regime II). 
In regime II, the shear rate is seen to decrease with stress for a controlled stress measurement, 
which supports the idea of coexistense between the SIS and the dilute phase at these shear 
rates §. Next, we observe a shear-thickened regime where the stress increases almost linearly 
with shear rate (Regime III). In this regime, the SIS forms a percolating cluster that undergoes 
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plug flow under shear and nucleates homogeneously in the sheared solution |3[] . Regime IV is 
a shear thinning regime arising out of the fracture of the SIS at very high shear rates. 

Figure 2 shows the shear stresses measured on application of step shear rates in the regimes 
II, III and IV. The shear stresses show an initial increase followed by a chaotic time-dependence 
about a mean value after a time t p . We would like to point out at this stage that the flow 
curve shown in Fig. 1 is not the true steady state flow curve and the data in regimes II and 
III are in all probability, transient. We come to this conclusion as we find from Fig. 2 that 
steady state is not reached even after shearing for 1500 seconds. We find that t p ~ j~ a , where 
a=1.8, in contrast to a=l for dilute TTAA/NaSal solutions ||. The emphasis of this work 
is on the analysis of the acquired time-series using standard non-linear dynamical methods. 
We prove the existence of deterministic dynamics in the stress relaxation of sheared CTAT for 
shear rates lying in regimes II and III, where we have derived finite values of the correlation 
dimensions and positive values of the Lyapunov exponents. The dynamics in regime IV shows 
a transition to higher degrees of complexity due to the stronger spatial dependence of the 
flow fields and a consequent coupling of the different wavevector modes, arising out of the 
percolation and fracture of the SIS. 

Data Analysis. We have employed the standard tools of nonlinear dynamics to 
calculate certain metric and dynamical invariants |^| characterising the time-series measured 
in our experiments and displayed in Fig. 2. This is done to distinguish whether our signals arise 
out of stochastic or deterministic processes. As opposed to stochastic noise, chaotic dynamics 
may be predicted over short time scales and are extremely sensitive to initial conditions. The 
first invariant we have calculated is the correlation dimension v which gives us an idea of the 
geometry of the attractor in phase space on which the stress trajectories lie in the asymptotic 
limit. Next, we have calculated the largest Lyapunov exponent A which provides a measure of 
the divergence of neighbouring trajectories. Presence of chaotic dynamics in a system requires 
v > 2 and A > 0. The analysis is done by reconstructing the phase space by embedding 
the experimental time-series of shear stress a in m dimensions using time delay vectors L |^| , 
so that the invariants associated with the dynamics may be calculated unambiguously. Let 
°~j — vQ) At) denote the time-series shown in Fig. 2, consisting of stresses measured at regular 
time intervals At (=ls), with j=l to N (N=1500), at various shear rates maintained constant 
during the duration of the experiment. A reasonable value of L ( =L ) has been estimated for 
each time-series by plotting the stress autocorrelation function C(r) versus r. The value of r 
at which C(r) falls to \ of its maximum value of 1 (at t — 0) is used as an estimate of the 
optimum delay time L Q . An optimal value of the minimal embedding dimension m Q has been 
estimated using the method of false nearest neighbours proposed by Kennel et al. |l(j, juj. 
The values of L thus obtained are used to construct the sequence of vectors {Xi : i=l,2,....N- 
(m-l)i } which describe trajectories in m-dimensional phase space arising out of the dynamics 
F : Xi — *• Xi+i, and which are used for the computation of the invariants characterising 
the dynamics such as the correlation dimension v and the Lyapunov exponent A ||. The 
calculations of the correlation dimension v uses the algorithm of Grassberger and Procaccia 
[ |L2"| which involves the calculation of the correlation integral C(R). The correlation integral 
C(R) is defined in an m-dimensional phase space as C(R) — J™^^ ^2 1=1 H(R — \X i —Xj\), 

where H(x) is the Heaviside function and \Xi — Xj | is the distance between the pair of points 
(i,j) in the m-dimensional embedding space. The sum in the expression for C{R) gives the 
number of point pairs separated by a distance less than R. For small Rs, C(R) is known to 
scale as C(R) ~ R" , where v is the correlation dimension A plateau in the plot of v versus 
log(i?) gives the correct v for a chosen embedding dimension m. If the attractor is unfolded by 
choosing a large enough m, then the correlation dimension v becomes independent of the value 
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of the embedding dimension to. The value of to = m for which this independence sets in is 
the correct embedding dimension, and the corresponding v is the correlation dimension. If an 
experimental signal satisfies v < to, then the signal is due to deterministic chaos rather than 
random noise |T^ ]. Figure 3 (a) shows a plot of v versus log(i?) for to=1-6 (increasing in the 
direction shown by the arrow) at 7=60s~ 1 , which corresponds to regime II of the flow curve. 
The value of is, obtained from the plateau in the v versus log(i?) plot, is found to saturate to a 
value of 2.6, which gives to ~ 3 or 4. The inset shows the values of the correlation dimensions 
v versus the embedding dimensions m calculated from the experimental data (open circles, 
marked as E) and five sets of surrogate data (SI-SV) generated by randomising the phases of 
the original data JlO| . The values of v calculated from the surrogate data do not converge even 
at to=6 and lie well above the values calculated from the experimental time-series. This implies 
the presence of deterministic dynamics in the stress relaxation of shear-thickening CTAT at 
7=60s _1 . 

Figure 3(b) shows the plot of v versus log(i2) for the data set acquired at 7=400s -1 which 
saturates to an approximate plateau value of 4.5 for m >6 over a range of log(i?). The values 
of is calculated from the experimental and surrogate time-series are not very different in this 
case. These large values of v and to calculated for 7 = 400s -1 may be interpreted as a 
manifestation of the increased complexity of the dynamics of the stress relaxation in regime 
IV arising out of the fracture of the SIS which is accompanied by the formation of vortices in 
the sheared solution. Figure 4 shows a plot of v versus 7 for the regimes II, III and IV. Figure 
5 shows the calculation of the largest Lyapunov exponent for 7=40, 50 and 90s -1 , following 
the algorithm of Gao and Zheng Q. This is done by defining dij(k) = \\X i+ k — Xj + k\\, 
the Euclidean distance between two vectors constructed using the embedding theorem by k 
iterations of the dynamics F, and plotting A=< ln[dij(k)/dij(0)] > as a function of k. The 
Lyapunov exponents, calculated by using the relation A = S'/Atln(2), where S is the slope 
of the plot, are found to be positive at all shear rates in regimes II and III, which confirms 
the existence of deterministic chaos at these shear rates. Due to the limited extent of our 
time-series, we have been unable to obtain reliable values of A for shear rates lying in regime 
IV. Preliminary analysis indicates an increase in the values of A in regime IV from those 
obtained in regime III, confirming the existence of higher dimensional chaos in this regime. 
The inset of Fig. 5 shows the values of A at different shear rates in regimes II and III. 

Conclusions. - In summary, we have proved the existence of deterministic chaos in the stress 
relaxation of dilute CTAT solutions in the shear-thickening regimes (Regimes II and III). The 
chaotic dynamics can arise out of the stick-slip between the SIS and the coexisting dilute phase 
in the sheared solution. In regime IV, where the spatial dependence of the flow fields becomes 
significant due to the percolation and fracture of the SIS, we have demonstrated an increase 
in the complexity of the dynamics, which is reflected by an increase in the optimal embedding 
dimension m , the correlation dimension v and the Lyapunov exponent A. Stress relaxation 
in dilute CTAT, therefore, provides yet another example of deterministic chaotic dynamics in 
a driven system arising out of the coupling between the lateral and normal degrees of freedom 
HQ. We hope that our work will stimulate theoretical understanding of the nonlinear flow of 
dilute surfactant solutions in terms of nonlinear coupled equations having additional predictive 
capabilities. 
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FIGURE CAPTIONS 

Fig. 1. - (a) shows a plot of the shear stress a while (b) shows a plot of the viscosity 77 
versus shear rate 7. Depending on the behaviour of 7/ at different shear rates, the flow curve 
may be divided into four regimes : regime I where the surfactant solution shows Newtonian 
flow behaviour, regime II characterised by shear-thickening as a result of the nucleation of 
SIS, regime III where the SIS spans the Couette and regime IV which shows shear-thinning 
as a result of the fracture of the SIS. The filled circles represent the shear-rate controlled 
measurements while the open squares denote stress control. 

Fig. 2. - Plots of the relaxation of the shear stress a at shear rates 7= (a) 40, (b) 50, (c) 
60 s- 1 (regime II), (d) 90, (e) 200 s" 1 (regime III) and (f) 400s" 1 (regime IV). 

Fig. 3. - The calculations of the correlation dimensions from the time-series of stresses 
at 7= (a) 60 _1 and (b) 400s _1 . The insets show the calculations of v for the experimental 
data set (denoted by E) and five sets of surrogate data ( SI-SV) generated according to the 
algorithm of Thcilcr et al. H]. v saturates to 2.6 for the experimental data at 7 = 60 _1 . At 
7 = 400 -1 , the z/s calculated from the experimental data and the surrogate data sets are not 
very different, which indicates a transition to higher-dimensional dynamics. 

Fig. 4. - Plot of v versus 7, which shows a monotonic increase of v versus 7. 

Fig. 5. - Calculations of A versus k for 7=40, 50 and 90s -1 . The Lyapunov exponents A 
which may be estimated from the slopes of the plots, are found to increase with 7, as shown 
in the inset. 
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